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We formulate a direct and systematic perturbation theory on the shift of the main paramagnetic 
peak in electron spin resonance, and derive a general expression up to the second order. It is applied 
to the one-dimensional XXZ and transverse Ising models in the high field limit, to obtain explicit 
results including the polarization dependence at arbitrary temperature. 



Introduction — Electron spin resonance (ESR) is not 
only an important experimental method, but also poses 
fundamental problems in theoretical physics. At the 
same time, it provides a good testing ground for the- 
ories, as very precise spectra can be obtained in many 
experiments. This is evident historically in the fact that 
the general formulation of linear response theory, the so- 
called Kubo formula, appeared first in a treatment of 
magnetic resonance by Kubo and Tomita, 1] and this be- 
came a standard theory on ESR. In addition to steady ad- 
vances in experimental techniques, a recent development 
of a field theory formulation 2, |j| for one-dimensional 
systems and also of numerical approaches ^ IH 1^ re- 
newed interest in ESR. 

In ESR, we study a magnetic system under an applied 
field H . The direction of the applied field is defined as z 
axis. The Hamiltonian is given by 



H — Hi 
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where 'Hsu{2) represents the isotropic interactions such 
as the Heisenberg exchange, and TL' is the anisotropy 
which is often small. TLz = —HS^ is the Zeeman term, 
and the notation 5" is used for the total spin operator 
as 3°" = Ylij Sj throughout this Letter. 

An electromagnetic radiation, with the frequency w of 
the same order as H, is applied and its absorption is 
measured. We consider the Faraday configuration, in 
which the oscillating magnetic field is perpendicular to 
the z axis. Generally, within the linear response the- 
ory, the absorption intensity is related to the imaginary 
part Xaa(^) of the complex dynamical spin susceptibility 
(i.e. the retarded Green's function of the spin operators) 
Xaai'-^) at zero momentum, where a is a spin axis deter- 
mined by the polarization of the electromagnetic wave. 
If the oscillating magnetic field is parallel to the x (y) di- 
rection, the absorption intensity is related to Xxx iXyy)- 
In the presence of general anisotropy, and S'^ are not 
equivalent, since they give different spectra. We call this 
dependence polarization dependence 0, 0, ■ 

In this Letter, for simplicity, we discuss x'aai^) ^ 
function of w at a fixed H, although usually in experi- 
ments the absorption intensity cx ujx'aai^) is measured 
with varying H while w is constant. When the spins are 
non-interacting, Ti = Tizi we trivially obtain 



We thus find the delta-function peak, which is called the 
paramagnetic resonance, at lo = H. A characteristic of 
ESR, which measures the zero wave vector component 
of x"i is that the lineshape still keeps the same form 
even in the presence of 'Hsu(2)- Nontrivial changes in 
the lineshape only occurs in the presence of anisotropy. 
For a small anisotropy 7i', the paramagnetic resonance 
aX uj = H \s often shifted and broadened. It is thus useful 
to obtain the shift and width of the main paramagnetic 
resonance peak in terms of the perturbation theory in 
the anisotropy Ti' , with the unperturbed Hamiltonian 

■^0 — T^z + 'Hsu{2)- 

However, this is not as easy as it may sound, even 
in a weakly interacting system. A realistic, broadened 
lineshape cannot be obtained in any finite order of the 
perturbative expansion of x" ■ This is apparent as x" h^^s 
the delta function at the zeroth order, while the "true" 
x" with a broadened resonance generally does not con- 
tain a delta function. The delta function could only be 
eliminated by a summation of the perturbation expansion 
of x" up to the infinite order. 

Several existing approaches to the shift and width are 
based on an infinite summation of the perturbation se- 
ries, relying on various grounds. In the Kubo- Tomita 
theory the infinite summation is carried out by as- 
suming that higher order terms in the cumulant expan- 
sion vanish. Although this assumption appears to be 
correct in certain cases, its range of validity has not 
been clarified. Indeed, it was pointed out recently that 
a straightforward application of the Kubo- Tomita theory 
fails in the presence of the Dzyaloshinskii-Moriya inter- 
action For S = 1/2 chain at low temperature, 
a field theory approach based on bosonization technique 
is developed 0)13 ■ There, the infinite summation is per- 
formed systematically in terms of the Feynman-Dyson 
self-energy. However, this is only possible because of the 
particular mapping of the ESR problem to the Green's 
function of the boson field, which holds only for one- 
dimensional system at low temperatures. Therefore, no 
immediate generalization of this technique to other cases 
of interest is possible. 

In the Mori-Kawasaki theory , an ansatz is made for 
the dynamical spin susceptibility as 
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in terms of the memory function M{oj), which is assumed 
to be smooth a,t uj ^ H. It is actually of the same form 
that results from the self-energy summation. With this 
ansatz, it is possible to determine the shift and width 
from the low order perturbative expansion of the suscep- 
tibility. However, it should be noted that in the absence 
of the diagrammatic perturbation theory, there is gener- 
ally no justification of the memory function ansatz lpj)l. 
Thus, its range of validity has not been clarified. 

The shift and width, on the other hand, have also been 
discussed in terms of the moments. The ESR shift may be 
related to the normalized first moment that is the "mean 
frequency" of the entire ESR absorption. However, it of- 
ten happens that, in the presence of the anisotropy, the 
peaks of the absorption spectrum appear at higher fre- 
quencies such as around 2H, 3H, AH .... In the presence 
of such satellite peaks, in order to discuss the shift of 
the main paramagnetic resonance to higher orders, one 
must calculate the projected moments which are defined 
by restricting the domain of integration to the support of 
the main paramagnetic resonance. However, the existing 
calculations [H O IH 111 III III 113 of the projected 
moments involve nontrivial assumptions, and their range 
of validity is again questionable. The polarization depen- 
dence has not been discussed in terms of these approaches 
based on the projected moments. 

Given this situation, 50 years after the publication of 
the Kubo-Tomita theory, it would be worthwhile to de- 
velop an alternative approach to ESR which is comple- 
mentary to the existing ones. 

Formulation — In this Letter, we aim to develop a per- 
turbation theory directly on the shift of the main para- 
magnetic resonance. 

For simplicity, we only consider a linear polarization 
of the oscillating magnetic field. When the oscillating 
magnetic field is parallel to the a axis, the ESR peak po- 
sition u!m,a may be defined by a maximum of x"^, which 
satisfies 



9ujXaai^m,a,a) = 0, 



(4) 



where a is the small perturbation parameter which gives 
the magnitude of the anisotropy H' . For the sake of 
brevity, we will omit the polarization index a from ujm^a 
until we discuss the polarization dependence. 

Our strategy is to expand the peak frequency in a as 



(5) 



However, an attempt to solve condition perturba- 
tively encounters a problem already at the zcroth order. 
This is because, at the zeroth order, the dynamical sus- 
ceptibility is given by the singular delta function as in (j^J. 
In fact, at the zeroth order, Q holds for any ^ H. 

The delta function can be defined as 



where e is often called the convergence factor, which is 
usually regarded as positive infinitesimal. The introduc- 
tion of the convergence factor corresponds to multiplying 
e~'^* with the retarded Green's functions, as functions of 
real time t. Generally, the inclusion of the convergence 
factor amounts to replacing uj with uj + ie. 

In order to circumvent the above problem due to the 
singularity, we perform the perturbative expansion 
keeping the convergence factor e finite. The limit e —>■ +0 
is taken after the expansion is obtained. In fact, a finite 
e is often used to estimate the continuous spectral func- 
tions from a finite-size numerical calculation, which in- 
evitably yields a collection of delta functions. Sometimes, 
e may be considered as a small but finite parameter that 
represents the relaxation effects which do not appear in 
the microscopic model but actually exist in a real system. 

If the shift is Taylor expandable in both e and a, this 
approach would give the correct perturbation series in 
a after taking the e ^ limit. This should be the 
case when the Mori-Kawasaki memory function ansatz 
is valid, as the shift is given by the real part of the mem- 
ory function, which is assumed to be smooth at w ~ 
and to be Taylor expandable in a. As we will see below, 
our formulation is actually more general than the Mori- 
Kawasaki ansatz. Namely, our approach gives correct 
answer even when the line shape is Gaussian. In con- 
trast, the ansatz Q does not match the Gaussian line 
shape. 

On the other hand, when the main paramagnetic reso- 
nance splits into multiple peaks, there would be multiple 
solutions to eq. at a; ~ _ff . Thus the present argument 
(as well as the Mori-Kawasaki memory function ansatz) 
does not apply to such cases in a straightforward manner. 
The generalization to this class of problems is discussed 
elsewhere this Letter, we focus on the simplest 

case with the single main paramagnetic resonance peak. 

Equation Q is then expanded as 

+a„9^£^(^(o)^ 0)^(1) + ldld^xla«, 0)} + • (7) 

Solving eq. {Tj) order by order with respect to a, we obtain 
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To evaluate these expressions, it is useful to derive the 
exact identities such as 



X+-(w) 



N{uj - H 



H + iey 



(10) 

where N is the total number of spins in the system, e is 
the convergence factor mentioned above, xh is the com- 
plex dynamical susceptibility of and S~ , and Xaa^ is 
that of A and A'' where A = [H', S+]. The identity ^ 
follows from the Heisenberg equation of motion 



dS+ 
dt 



iA, 



(11) 



and its Hermitian conjugate (see also Appendix of 
ref. 01 )• As discussed above, we regard the convergence 
factor e as being finite until the expansion ||SJ) of the shift 
is obtained. 

After straightforward but tedious calculations, we ob- 
tain the shift of the main paramagnetic resonance up to 
the second order. The first order result reads 



{[AS-])o 

2(5-)o ' 
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where we have introduced the notation ()„ as the n-th 
order contribution (in a) to the expectation value defined 
with respect to the full Hamiltonian This result is 
independent of the polarization. Although it is found 
to be identical to the well-known formula by Kanamori 
and Tachiki and others [l^I , our derivation is more 
systematic than the previous ones. 

The second order result, on the other hand, is more 
involved and interesting as shown below. Unlike in the 
first order term H12|l. the shift in second order generally 
has a polarization dependence. Thus here we restore the 
polarization index a in the shift 5ujm,a, which is given as 

{[A,S-])o + {[A,S-]h 



— — — lim 



(S^), ^ {[A^,S-])o-{[A,S+])o 



± 



\ (e [xaa + Xa-<a^) - {Qu^x'aa + '9^x!4Ut)) 

(13) 
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where x' is ttiG real part of the complex dynamical sus- 
ceptibility, x' Slid x" evaluated at uj — LUm^ , and the 
double sign ± takes -I- or — for a = x and y, respectively; 
a is the polarization direction of the oscillating magnetic 
field. The term "0{e^)" refers to the contributions with 
prefactors explicitly containing or higher powers of e. 



Because of the limit e +0, the terms explicitly multi- 
plied by powers of e, including the "0(e^)" terms, appear 
to vanish. However, the correlation functions in eq. H13|) 
can have singularities as powers of 1/e, leaving a finite 
contribution together with the explicit powers of e. Thus 
we cannot ignore such terms a priori, and we show later 
a concrete example with such a nontrivial contribution. 
For brevity, in eq. H13(l we have omitted writing out the 
"0(e'^)" terms explicitly, since they actually vanish in all 
the examples we have studied so far. 

Eq. ifT^ should be valid up to the second order in the 
perturbation parameter a, in any dimension, at any tem- 
perature, for any spin, and for any interaction strength. 
However, it must be noted that the evaluation of the cor- 
relation functions appearing in cq. H13|l is still a nontrivial 
problem. 

Applications — As a demonstration, here we apply our 
general second-order result (|13() to the high field limit 
Hz ^ 'H.su{2),'H' . In this case, the isotropic interac- 
tion 'Hsu(2) (if any) may be included in the perturbation 
H' . This is convenient for calculation, because the un- 
perturbed Hamiltonian Tip = Tiz is trivially solvable. 
As a consequence, eq. (fT^ can be evaluated exactly. Al- 
though this is a weak-coupling limit, the problem of the 
ESR lineshape is by no means trivial, as is evident in the 
results to be shown. Here we study a few one-dimensional 
models, for which we can verify our results by comparing 
them to known exact solutions in certain limits. 

First we consider the 5 = 1/2 XXZ chain in the high 
field limit. Here, as discussed above, all the interactions 
including the isotropic part are included in the perturba- 
tion: 



i+i' 



(14) 



where J,J'^H and z is assumed to be identical to the 
direction of the applied magnetic field. For this model, 
in eq. (|13|l the terms involving dynamical susceptibility 
turn out to vanish, and contributions from the static ex- 
pectation values give the final result: 



j'tanh 



J 



,2 /3coth(/3i7) 



2 J ' " 1 -I- cosh (/3F) 
-JJ'^sech^r^Vanhr^ 



V 2 



,(15) 



where /3 is the inverse temperature. Because of the U(l) 
rotational symmetry about the z axis, there is no polar- 
ization dependence in this model. 

In the zero temperature limit ^ oo, eq. (|15|l reduces 
to J'. In this limit, the ESR probes an excitation gener- 
ated out of the ground state. Since we consider the strong 
field case H ^ J, J', the ground state is simply the fully 
polarized state and the elementary excitation is given by 
a magnon corresponding to single spin flip. Thus, in the 
zero temperature limit, the ESR spectrum consists of a 
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single delta function at oj ~ H + J' , corresponding to 
the magnon energy at the zero wave vector. Namely, the 
shift is J', exactly as obtained in our approach. 

In the infinite temperature limit /? — > 0, eq. I|15|) re- 
duces to J' /{2H). We emphasize that the ESR line- 
shape in the infinite temperature limit /3 — > is a non- 
trivial problem, although all the static correlation func- 
tions become trivial. In fact, to our knowledge, no exact 
result is available for the general J/ J'. Nevertheless, for 
the special case of J = J', the model reduces to the 
XX model, in which the infinite temperature limit of the 
lineshape is known exactly [Ull^. (See also ref. for 
comparison to the Kubo-Tomita theory and to a numer- 
ical calculation.) We confirmed that the shift obtained 
from this exact lineshape agrees with our result above. 
Thus, our result l)15|l is verified at both the zero and 
infinite temperature limits. We believe that eq. 1)15(1 is 
exact as the second order perturbation expansion at any 
temperature. 

Next we turn to the transverse Ising model in the high 
field limit. We regard the Ising interaction 



nj/ t/ \ ^ QX QX 



(16) 



as the perturbation of the Zeeman term. Note that the 
quantization axis of the Ising interaction is the x-axis and 
thus the model is the transverse Ising model. This model 
lacks U(l) rotational symmetry about the z-axis. Thus 
it is possible to have a polarization dependence, which 
we have indeed obtained. Technically, the evaluation of 
eq. 1(131) for this model is more involved than in the XXZ 
model discussed above, as several dynamical correlation 
functions such as 'x!aa non-vanishing contributions 
to the shift. In particular, we find nonvanishing contri- 
butions from the terms 

e {xaa + Xa^a^ ) - {du^XAA + d^XA^A^ ) (17) 

in the e — > limit, due to the singularity in the dynamical 
correlations. Collecting all of the contributions, eq. 113(1 
for the transverse Ising model reads 



5uJn 



= — tanh 



f) 



'1 -f- tanh 



f) 



A(3H 



al3H 



-I3H 



(18) 



J'2 

'3 — 3 tanh' 



if 
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P0H 



-I3H 



(19) 



for the linear polarization parallel to the x and y axes, 
respectively. The nontrivial polarization dependence is 
evident in the result. 



In the infinite temperature limit /3 ^ 0, they are fur- 
ther simplified as 5oJm,x J'^/(16i?) and Sujm.y — > 
5J'^/{16H). Namely, there is a difference of factor 5 
between the different polarizations. They, including the 
polarization dependence, are again found to agree with 
the known exact result |2l| on the transverse Ising model 
in the infinite temperature limit. In this case, the polar- 
ization dependence, which is necessary for the agreement 
with the exact result, comes entirely from the nontrivial 
contribution fT7|l . 

On the other hand, our exact second-order results are 
obtained for the entire temperature range. In contrast, 
the Kubo-Tomita theory and its modification assumes 
a high temperature, and the previous approach [13 . Iisl 
01 based on the projected moment fails to give the po- 
larization dependence. 

Summary — To summarize, we have developed a direct 
perturbation theory on the shift of the main paramag- 
netic resonance in ESR, and have obtained a general and 
explicit formula up to the second order. As a demon- 
stration, it is applied to the S — 1/2 XXZ chain and the 
transverse Ising chain in the strong field regime, and the 
shift in the second order is obtained in closed forms for 
the entire temperature range. The agreement with the 
exact results available in a few limiting cases is verified. 

Our method is also applicable to systems with higher 
spins and in higher dimensions. The evaluation of the 
formulae in these cases is possible, at least in the high- 
field limit. Furthermore, even for strongly interacting 
systems, we could utilize various methods such as field 
theory and exact solutions to evaluate eq. ((13(1 . 

We also expect that the present approach can be ex- 
tended to the problem of the ESR width, as well as other 
problems of interest in quantum dynamics. More details 
and further generalization will be given in a separate pub- 
lication llTl. 
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